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Abstract. We report on the collapse and revival of Ramsey fringe visibility when a spatially dependent 
phase is imprinted in the coherences of a trapped ensemble of two-level atoms. The phase is imprinted via 
the light shift from a Gaussian laser beam which couples the dynamics of internal and external degrees of 
freedom for the atoms in an echo spectroscopy sequence. The observed revivals are directly linked to the 
oscillatory motion of atoms in the trap. An understanding of the effect is important for quantum state 
engineering of trapped atoms. 
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1 Introduction 

A trapped gas of atoms can act as dispersive medium 
with a refractive index which depends on the internal 
state of the atoms. The state-dependent phase shift of 
probe laser light propagating through an ensemble of Cs 
atoms has recently been used to observe Rabi flopping 
on the clock transition non-destructively |l|2j . Such non- 
destructive measurements of a collective atomic quantum 
state component holds the promise to predict the outcome 
of subsequent measurements beyond the standard quan- 
tum limit [3^. This reduction in uncertainty is referred 
to as conditional squeezing and the resulting nonclassi- 
cal atomic state may be used to increase the precision of 
atomic clocks ^ 

In a recent paper [5] , we considered the effect of inho- 
mogeneous light shifts on the atomic quantum state evolu- 
tion when using a Gaussian laser beam for dispersive prob- 
ing of an ensemble of atoms confined in a dipole trap. The 
spatial intensity distribution of the probe beam implies 
that an atom experiences a position dependent differen- 
tial ac Stark shift of the clock levels and the atomic cloud 
acquires a spatial phase imprint as illustrated schemat- 
ically in Fig. [T] In the present paper we shall focus on 
the fact that the individual atoms are not stationary in 
the trap, but move about to explore regions of different 
probe light intensities. This is of consequence for how well 
their inhomogeneous phase spread can be compensated 
for by using Hahn echo techniques. Specifically, we shall 
investigate the degradation of Ramsey fringe contrast in 
echo spectroscopy as a result of atomic movement in be- 



Atoms 




Spatial Phase Imprint 



e-mail: kj aergaardQnbi . dk 



Gaussian Laser Beam 



Fig. 1. The Gaussian intensity distribution of a laser beam is 
mapped onto the atoms as a spatially dependent phase of the 
clock state superpositions. 



tween two perturbing light pulses on either side of the 
echo pulse. The fringe contrast can be observed to revive 
at half-integer multipla of the radial trap period and we 
demonstrate how this effect can be used to measure the 
trap frequency in situ without actually exciting collective 
oscillation modes. The specific form of the collapse and 
revival of fringe visibility is modelled readily when the 
anharmonicity of the trapping potential is taken into ac- 
count. 



2 Ramsey Interrogation and Echo 
Spectroscopy 

2.1 Ramsey spectroscopy of two-level atoms 

In Ramsey spectroscopy, a collection of two-level atoms 
with quantum levels 1 1) and 1 1) separated by an energy 
difference hQo interacts with two near resonant fields of 



Daniel Oblak et al.: Echo Spectroscopy of Atomic Dynamics in a Gaussian Trap via Phase Imprints 



nil 



light 



light 



Till 



— MW- 



E 



1 



ti 



-►-4- 



t. 



Fig. 2. Echo spectroscopy time hne as employed in our ex- 
periments. A near-resonant microwave (MW) field is applied 
to the atoms in a 7r/2 — tt — 7r/2 pulse sequence. The effect of 
short perturbing light pulses occurring at a time ti before and 
a time t2 after the vr-pulse is investigated. 



Several studies focus on the influence of the inhomoge- 
neous dipole trapping light and it has been demonstrated 
that an echo pulse is efficient in compensating the effect of 
the trapping laser on a time scale <C h/AE^ where AE is 
the differential ac Stark shift between | [) and 1 1) [10 . Fur- 
thermore, it has been shown that the differential light shift 
for an oscillator mode of a far-off-resonant trap (FORT) is 
well described by its time averaged value [9 such that in- 
homogeneous dephasing introduced by the Gaussian trap- 
ping laser beam profile can in many cases be reversed in 
an echo sequence. 



angular frequency Q ^ Qq separated by a time T ^6j . With 
all the atoms initially in the state | [) the first field inter- 
action — a 7r/2-pulse — produces an equal coherent su- 
perposition state (I i) + I ]))/\/2 for each atom. The atoms 
now evolve freely for the time T during which a phase (j) = 
(i7-i7o)T is accumulated by the state: (I i) + e''t'\]))/^/2. 
Finally, a second 7r/2-pulse is applied and the population 
difference of the two states is measured. This quantity will 
vary periodically with (j) giving rise to so-called Ramsey 
fringes when either T or i? are scanned. 

If the two energy levels | [) and | j) are perturbed dif- 
ferentially in any way during the free evolution period, 
the energy splitting HQq and thus (j) will be affected. This 
may occur homogeneously such that all the atoms are per- 
turbed by the same amount causing an over all shift of 
the Ramsey fringe pattern [7] . Essentially, this fringe shift 
constitutes an interferometric measurement of the pertur- 
bation strength. For nonuniform perturbations the phase 
(j) differs from atom to atom and in the Ramsey measure- 
ment of the whole collection of atoms the fringe shift is 
accompanied by a degraded visibility as a result of inho- 
mogeneous dephasing [5]. 



2.2 Coherence echoes 

Effects of any form of inhomogeneous dephasing encoun- 
tered in Ramsey spectroscopy can, to a large extent, be 
cancelled by introducing a so-called echo-pulse (yr-pulse) 
between the two 7r/2-pulses at time ti in the Ramsey se- 
quence [8 . The echo pulse essentially inverts the sign of 
the phase (j) accumulated so far by each atom and if sub- 
sequently each atom encounters the same amount of per- 
turbation as prior to the 7r-pulse a rephasing will occur at 
time ^2 = h after the echo pulse. In this sense the per- 
turbation is reversible e.g. for a collection of stationary 
atoms in a inhomogeneous light field. Irreversible dephas- 
ing may result from fluctuating perturbations such as the 
variations in phase shift due to noise in the intensity of the 
inhomogeneous light field or movement of atoms therein 

i- 

Since echo spectroscopy nullifies the effect of reversible 
perturbations it has proven to be a powerful tool for inves- 
tigating the nature and magnitude of atomic decoherence 
due to the combined effect of irreversible processes includ- 
ing dephasing and spontaneous decay of atoms |9|1Q|11|12| . 



2.3 Inhomogeneous light shifts and trap dynamics 

In the present paper we shall consider the perturbations 
from an auxiliary inhomogeneous pulse of light applied 
before the echo pulse and an identical light pulse applied 
after the echo pulse as indicated in Fig. [2] For a stationary 
atomic ensemble perfect rephasing would be expected ex- 
cept for decoherence of the atomic state due to irreversible 
spontaneous scattering events and quantum fluctuations 
(shot noise) of the two light pulses. Hence, a measurement 
of the Ramsey fringe visibility would appear to be ideally 
suited for determining the amount of decoherence intro- 
duced by a given laser beam. The situation is somewhat 
complicated when atoms change positions between the ap- 
plication of the two light pulses. To elaborate on this issue 
we consider a simple ID model of particles evolving in a 
harmonic trap. The distributions of atoms in momentum 
and poisition are represented as a dynamical phase space 
plot [see Fig. [s] (a)]. The application of a light pulse with 
Gaussian intensity profile will imprint a phase which only 
depends on the position if we assume that the duration of 
the light pulse is much shorter than an oscillation period 
T in the trap. Now, during the free evolution the imprint 
will rotate in dynamical phase space at angular frequency 
CO = 2it/T. In Fig. IsFb), we outline the situation for echo 
spectroscopy in the cases (i) n = r2 = 7r/4a; and (ii) 
Ti = T2 = 7r/2cc;, where n and T2 are the durations from 
the light pulse before and after to the echo pulse, respec- 
tively. Complete rephasing of the atoms is encountered 
only if ri+r2 = 7r(mod 7r)/cc;, i.e. the atoms oscillate in the 
trap for a half-integer multipla of a trap period between 
the two light pulses. For an echo spectroscopy experiment 
this implies that the Ramsey fringe contrast will depend 
on the time separation ri + T2 between the light pulses. 



3 Experimental 

3.1 Setup 

Details of our experimental setup and atomic sample prepa- 
ration can found in ^ . The starting point for the exper- 
iment presented in the present paper is an ensemble of 
- 10,000 - 50,000 Cs atoms polarized in the 65'i/2(F = 
3,7Tii? = 0) = ID clock state and confined by a ~ 4 W 
Yb:YAG laser beam focussed to a waist of ~ 40/im. This 
dipole trap is characterized by oscillation frequencies in 
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Fig. 3. (a) Dynamical phase space representation of atoms in a ID harmonic oscillator potential of angular frequency uj. (b) 
Phase imprint on the atoms from the position dependent light intensity of Gaussian laser beam during echo spectroscopy, where 
two identical light pulses act on each side of the echo vr-pulse. (i) In general echo spectroscopy cannot completely rephase the 
atoms, (ii) In the special case of the time between the two light pulses being n -\- T2 — 7r(mod7r)/a;, the effects of the two light 
pulses cancel. 



the order of a kHz radially and a few Hz axially. The 
clock transition 1 1) ^^ 651/2 (F = 4, m^ = 0) = | j) is 
driven using microwave radiation around 9.2 GHz applied 
in a 7r/2 — tt — 7r/2 echo spectroscopy sequence as depicted 
in Fig. |2j At the end of this sequence we determine the 
fraction of atoms residing in the 1 1) state. The probing of 
atoms is performed with a beam of light propagating along 
the trap axis with a waist of 18 jiia located at the center 
of the trap. The frequency of the probe light is blue de- 
tuned by 160 MHz from the 66'i/2(F = 4) ^ 6P3/2(F = 5) 
transition and via the dispersive atom-light interaction the 
probe light experiences a phase shift proportional to the 
number of 1 1) -atoms which is measured using a shot noise 
limited Mach Zehnder interferometer [4,13j . By applying 
light resonant with the 66'i/2(F = 3) ^ 6P3/2(F = 4) 
transition all atoms are pumped into the 65'i/2(P = 4) 
level and the total number of atoms involved is determined 
from a subsequent phase shift measurement. Optionally, 
we can apply light from an additional probe laser which 
is red detuned by 135 MHz from the 6S'i/2(F = 3) ^ 
6P3/2(P = 2) transition and hence couples to the | [) pop- 
ulation. By engaging the two probes simultaneously we 
can obtain a zero (mean) interferometer phase shift for 
ensembles in an equal clock state superposition irrespec- 
tive of the total number of atoms. This two probe color 
configuration has proven convenient in our measurements 
of atomic projection noise ^. The two probes of this two- 
color scheme are merged in a single mode optical fiber to 
ensure good spatial overlap and hence enter the interfer- 
ometer at the same input port. 



tion gained and coherence lost. Establishing an optimal 
balance between decoherence and measurement strength 
is important for quantum state engineering of a squeezed 
clock state population difference via a dispersive measure- 
ment. As it turns out, the optical depth of the atomic 
sample is the key figure of merit determining the optimal 
optical power and thus the amount of decoherence [14]. 
Preferably, and to achieve a high degree of squeezing, the 
optical depth should be large such that each probe photon 
interacts with many atoms. In this respect the ~ 1 : 200 
radial to axial aspect ratio of our sample provided by the 
dipole trap potential is favorable and gives rise to an op- 
tical depth of up to ~ 20 along the direction of the probe 
laser beam. 

From these considerations, it is obviously important to 
have a handle on the amount of decoherence introduced 
by the dispersive probing scheme and echo spectroscopy 
would appear to be the method of choice [TT]. Ideally, a 
light pulse derived from the probe laser could be devided 
into two with each part applied before and after the echo 
pulse, respectively, and the reduction in Ramsey fringe 
amplitude would then gauge the decoherence from probe 
light considered as a perturbation. This method relies on 
the complete cancellation of the reversible dephasing by 
the inhomogeneous light shift. However, the movement of 
atoms in between the two pertubing pulses may also lead 
to imperfect rephasing causing a Ramsey fringe reduction 
as outlined in |2.3| In fact, this effect is significant and is 
prominently manifested in the echo spectroscopy. 



3.2 Decoherence and inhomogeneous dephasing 

For low optical powers the dispersive probing scheme is 
close to nondestructive in the sense that spontaneous scat- 
tering is very limited. However, the ratio of signal to noise 
(which in our case is the shot noise of light) for a phase 
shift measurements increases with increasing probe pho- 
ton number. Hence, there is a trade off between informa- 



3.3 Results 

3.3.1 Echo Spectroscopy 



As described in section 2.3 the dynamics of atoms in our 
dipole trap is expected to affect the ability of an echo 
pulse to balance the effect of two surrounding phase shift- 
ing light pulses, as can indeed be observed experimen- 
tally in echo spectroscopy. To illuminate the effect, we 
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present in Fig. |4] examples of Ramsey fringes as recorded 
in the two-color probing scheme using the echo sequence 
of Fig. [2] for light pulse separations ri + r2 of 12, 250, 
and 500 /is, respectively. The light pulses have a dura- 
tion of 2 jiis which is short compared to all other relevant 
time scales in the experiment such as the trap oscillation 
period of ~ 1000 jas. For the traces shown, ti was kept 
fixed at 1500 ms while ^2 was scanned between 500 ms 
and 2500 ms. To enable the observation of Ramsey fringes 
in the time domain the frequency of our microwave source 
was detuned by 3 kHz. The detuning sets the undulation 
frequency of the Ramsey fringes which have their maxi- 
mum amplitude at t2 = ti = 1500 ms due to the echo 
rephasing pulse. 

The most intriguing aspect of Fig. [4] is the collapse and 
revival of fringe visibility. For a time separation ri + r2 = 
12 /is between the light pulses which is short compared to 
the trap frequency, the echo pulse is effective in restoring 
the Ramsey fringe visibility: Since the light pulses are so 
closely spaced that the atoms hardly have time to move in 
between, the second pulse essentially undoes the inhomo- 
geneous phase imprint of the first. An upper bound on the 
amount of decoherence can be estimated from reduction in 
fringe amplitude, which in this case is 76%. A completely 
different situation is encountered at n + r2 = 250 /is with 
a light pulse separation in the vicinity of a quarter of a 
radial trap period. Here the atomic ensemble develops a 
nonuniform phase as indicated in Fig.[3fc)(i) and the Ram- 
sey fringes for individual atoms will in general not interfere 
constructively. Hence a certain degradation in the ensem- 
ble fringe contrast happens. Finally, for a pulse separation 
of Ti + r2 = 500 /is in the vicinity of half a radial trap pe- 
riod a Ramsey fringe revival happens. At the time of the 
second light pulse an atom will find itself close to the ra- 
dial distance (from the trap center) it had when the first 
pulse was applied and hence experience the same light in- 
tensity. Due to the inverting effect of the interposed echo 
pulse a phase shift cancelation will occur. 
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Fig. 4. Ramsey fringes as recorded in echo spectroscopy as 
a function of the time separation n + T2 between two light 
shifting pulses around the echo pulse. The fringe amplitudes 
have been normalized to the case when no light is applied. For 
increasing separation times the fringes are observed to first 
wash out and subsequently revive. 
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Fig. 5. Observed revival frequency of Ramsey fringes in echo 
spectroscopy as a function of the dipole trapping power. A 
square root scaling (line) describes the data well. Each data 
point was extracted from a revival trace as shown in the inset 
as recorded using the single color probing scheme. 



3.3.2 Revival frequency versus trap power 

The revival of Ramsey fringe visibility, resulting from in- 
terference between spatially imprinted phases in the atomic 
coherences, is linked to the radial trap oscillation frequency. 
Measurements of trap frequencies in cold atom experi- 
ments are typically performed by exciting either monopole 
or dipole oscillation modes, or by driving parametric losses 
when modulating the trap jl5|16j . Using the fringe revival, 
we are able to extract the trap frequency in situ without 
actually exciting motion — an atom is tagged in its inter- 
nal degrees of freedom according to its position. Towards 
this end we simply measure the height of the central Ram- 
sey fringe in our echo sequence Fig. [2] as a function of the 
separation between the two light shifting pulses. Hence, we 
keep ^2 = ti fixed and vary ri + r2. As the height of the 
central Ramsey fringe provides measure of the fringe vis- 
ibility this quantity oscillates in ri + r2 at twice the trap 
frequency as discussed above. Figure [5] shows the mea- 
sured revival frequency as function of the optical power of 
our dipole trapping beam. The observed revival frequen- 
cies are described well by a square root dependency on the 
optical power as expected from theory [15]. 



3.3.3 IVIodeling the revivals 

In a perfect harmonic trap the radial oscillation period 
is independent of oscillation amplitude and complete re- 
vival of the Ramsey fringe would hence be expected when 
Ti + r2 equals half a trap period. This is clearly not the 
case for the revival presented in the inset of Fig. [5] A pos- 
sible explanation could be sought in the anharmonicity 
of the Gaussian trapping potential. To investigate this 
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in more detail we perform numerical simulations of the 
phase evolution of N = 50, 000 particles in a Gaussian 
trap when subjected to the echo spectroscopy sequence 
shown in Fig. |2] including a non-uniform light shift from a 
Gaussian laser beam. During half a radial trapping period 
the axial movement of even the most energetic, trapped 
particles will only amount to a small fraction of the probe 
beam Rayleigh length which characterizes the length scale 
for the light-atom interaction volume. Furthermore, colli- 
sions between particles are expected to be negligible for 
the low densities and short time scales involved [17]. In 
addition, we shall assume that the mean kinetic energy 
for trapped particles is sufficiently large for the radial fre- 
quency variation along the trap axis to be negligible com- 
pared with the local frequency spread caused by the trap 
potential anharmonicity. Hence, as an approximation, we 
shall restrict our treatment to the 2D case of radial dy- 
namics. The details of the implementation of our numeri- 
cal simulations is given in appendix |AJ An expected fringe 
amplitude (relative to the unperturbed case) as function of 
time separation between the two perturbing light pulses, 
t = Ti + r2 , is given by 



nt) = 



Y^^_ g-p,'(*/2+t2)/7'»o, 



^Tit) 



E 



N 
j=0 < 



-p2(t/2+t2)/72^g 



(1) 



echo spectroscopy experiments for various perturbing light 
powers Ppert- We model our data using J^{t) as provided 
by our numerical simulations. The functional behavior of 
J^{t) depends on the three parameters T, 7, and 0o, i-e., 
the cloud temperature prior to the separatrix truncation 
(see appendix [A]) , the probe to trap beam ratio, and the 
peak phase shift. These parameters are adjusted to mini- 
mize the mean square difference to our experimental data. 
Least squares are obtained for kT/Uo = 1, 7 = 0.6, and 
0o/^pert = 5.0 X 10~^ which corresponds well to what 
would be expected in the experiment. Fixing the free pa- 
rameters of our simulation code to these values, we plot 
the normalized echo amplitude in Fig. [6] We note that 
our simple 2D model gives an excellent over-all agreement 
between experiment and simulation. 

The model does not include decoherence with the re- 
sult that all the simulated curves in Fig. |6] indicate perfect 
rephasing at t = 0. In the experimental data it is evident 
that rephasing is not perfect especially for the large probe 
powers. However, the effect of trap dynamics on the fringe 
visibility clearly dominates over the decoherence for most 
pulse separations and it is unclear how much the, in prin- 
ciple, reversible dephasing contributes to the fringe reduc- 
tion at the smallest pulse separation of 12/i5. Thus, mea- 
surements provide an upper bound on the (irreversible) 
decoherence 



(i) 

where pj and ^V^ are the radial position and echo phase 
[cf. Eq. Q], respectively, of particle j in the ensemble. 
Since the atomic ensemble is probed with a Gaussian light 
beam, each particle contributes to the total signal by a 
weight according to the intensity at its radial position 
at the time of probing. Figure |6] shows the fringe am- 
plitude versus light pulse separation time as recorded in 
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Fig. 6. Collapse and revival of the fringe amplitude in echo 
spectroscopy. Filled circles shows the experimental normalized 
fringe amplitude as function of the time separation n + T2 
between peturbing light pulses around the echo pulse for ten 
different values of the number of perturbing photons A^pert- 
Full lines show the result of our numerical simulations of the 
dynamics. 



4 Discussion 

In summary, we have reported on the manifestation of 
motional dynamics of trapped atoms in echo spectroscopy 
on the Cs clock transition. The movement of particles im- 
plies that the ac Stark shifts from two identical light pulses 
generally do not cancel completely on the application of 
an in-between Hahn echo. Rather, the echo fringe ampli- 
tude is observed to decrease with increasing separation of 
the two light pulses. For obtaining squeezing on the clock 
transition via dispersive measurements this is unfortunate 
for two reasons. First, the effect limits the applicability 
of echo techniques to restore the adverse inhomogeneous 
light shift effects of a strong off-resonant probe pulse. Sec- 
ond, using echo spectroscopy to gauge the atomic decoher- 
ence resulting from spontaneously scattered probe pho- 
tons when probing an ensemble is not straight forward. 
Obviously, an understanding of the effect of motional dy- 
namics is important when dispersive measurements are 
used for quantum state engineering. 

An alternative approach to eliminate inhomogeneous 
differential light shift of the clock states induced by a 
probe pulse is to apply a simultaneous light pulse from 
a second laser [18 . The frequency of the second probe 
laser should be chosen so that the differential ac Stark 
shifts introduced by each probe laser exactly cancel. In 
the present configuration, where the two probe beams en- 
ter the Mach Zehnder interferometer via the same input 
port, the detunings are pegged such that simultaneous ap- 
plication of the two probes on a coherent superposition 
state yields a zero mean phase shift. At these detunings 
the ac Stark shift imposed on each clock state add in sign 
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and instead of cancelling the differential light shift, the 
two-color probing scheme essentially doubles it. However, 
injecting the two beams via different ports of the inter- 
ferometer input beam splitter establishes a scheme where 
the detunings can be chosen to eliminate the differential 
light shift and achieve a state sensitive, balanced interfer- 
ometer signal at the same time [19 . We are presently in 
the course of reconfiguring our experimental setup to this 
favourable scheme. Echo spectroscopy experiments along 
the lines the present paper should then easily indicate if 
light shift cancelation has been accomplished. 

We recall that the primary motivation for our exper- 
imental efforts is to demonstrate squeezing on the clock 
transition. Here a balanced interferometer measurement of 
an ensemble in an equal clock state superposition will be 
used to predict the outcome of a subsequent measurement 
beyond the standard quantum limit |2|3|4j . The present 
paper and ref. [5 have investigated some adverse effects on 
the collective ensemble state which may accompany these 
"nondestructive" measurements and discussed routes to 
minimize them. Even when applying these strategies to 
eliminate inhomogeneous light shifts, the atomic trap dy- 
namics will lead to retardation effects when comparing 
two consecutive measurements on an ensemble. Particle 
motion and a Gaussian probe beam profile implies that 
an atom in general contributes to the interferometer sig- 
nal by different weights in the two measurements. A fuller 
treatment of this is beyond the scope of this paper, but it 
is expected that the degree of squeezing (i.e. the amount 
of correlation between the two pulses) is going the dimin- 
ish and increase analogous to the collapse and revival of 
echo fringe amplitude as reported on here. 
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A Methods for numerical simulations 

We obtain the initial conditions for our numerical simu- 
lations by considering a canonical ensemble of particles 
with mass tti at a temperature T. Using polar coordinates 
(p, ^), the phase space density is 

aip, 0,pp,pe) oc e-H+pl/p'y^^^+^^PW^T^ (2) 

where k is the Boltzmann constant and Pp and Pp are 
the conjugate momenta to p and 6>, respectively [20 . For 
the Gaussian trap the confinement potential is V{p) = 
—Uoe~^'^^ /^o)^ which has the (finite) trap depth Uq and 
depends only on the radial distance p to the trap center 
at the origin. We want to restrict ourselves to particles 
inside the separatrix of stable motion, i.e., we disregard 
unbound particles for which (p^ -^ p'g/ p^)/2m + V{p) > 
(see Fig.[7|. Integrating over particles inside the separatrix 




Fig. 7. Illustration of the Gaussian potential separatrix of 
bound motion projected onto {p,Pp,pe) phase-space. Particles 
inside the separatrix are trapped and form the starting point 
for our numerical simulations. 



yields a spatial density distribution 
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-Vip)/kT 
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exp I 



exp I 



.fZo V 



.2p! 



/7o>fcT 
Uo^kT 



(3) 

which for all practical purposes can be assumed to vanish 
for p ^ wq. Applying a finite cut-off radius po ^ u)o^ we 
assign random radii to the particles of our ensemble ac- 
cording to a probability density ex pe~^^P^^^^ using the 
rejection method [21 on the interval [0...po]- Next we as- 
sign the canonical momenta Pp and pe for each particle in 
correspondence with normal distributions of widths \/kT 
and pVkT^ respectively. We discard untrapped particles 
as described above and tag each trapped particle j with a 
phase (l)\^^ = ^oexp(— 2p^/7^k;q), where 7 is the ratio be- 
tween the waists of the trapping laser and the perturbing 
laser and (j)o is the peak phase shift. We finally integrate 
the equations of motion using the Verlet method to obtain 
pj (t) from which the echo phase 



4^-)^ 



'PTit) = '/'F^ - </'oexp(-2p2(t)/^2^2) 



(4) 



can be calculated. 
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